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If the solution to the solar neutrino puzzle falls in the LMA region, KamLAND should be able 
to measure with good precision the corresponding oscillation parameters after a few years of data 
taking. Assuming a positive signal, we study their expected sensitivity to the solar parameters 
(Si2,Am| 1 ) when considered in the framework of three-neutrino mixing after taking into account 
our ignorance on the mixing angle 6*13 . We find a simple "scaling" dependence of the reconstructed 
9i2 range with the value of #13 while the AmJ 2 range is practically unaffected. Our results show that 
the net effect is approximately equivalent to an uncertainty on the overall neutrino flux normalization 
of up to ~ 10%. 



The Sudbury Neutrino Observatory (SNO) measure- 
ment on the charged current reaction for solar neutrino 
absorption in deuterium jjj has provided an important 
piece of information in the path to solve the solar neu- 
trino problem (SNP) M JJ. In particular, all the post- 
SNO global analysis ]7[-pdJ have shown that the inclu- 
sion of the SNO results have further strengthen the case 
for solar neutrino oscillations with large mixing angles 
with best fit in the region with larger Am 2 , LMA. Un- 
fortunately, the LMA region is broad in mixing and mass 
splitting due to the uncertainties in the solar fluxes and 
the lack of detailed data in the low energy range. 

This situation will be improved in the very near fu- 
ture, in particular, if LMA is the right solution to the 
SNP. If Nature has arranged things favourably, the SNO 
measurements of the day-night asymmetry and the neu- 
tral to charged current ratio could help to identify and 
constrain the LMA solution jpU| . Furthermore, the ter- 
restrial experiment, KamLAND |12]| , should be able to 
identify the oscillation signal and significantly constrain 
the region of parameters for the LMA solution. 

The KamLAND reactor neutrino experiment, which is 
to start taking data very soon, is sensitive to the LMA 
region of the solar neutrino parameter space. After a few 
years of data taking, it should be capable of either ex- 
cluding the entire LMA region or, not only establishing 
Ve ^othor oscillations, but also measuring the oscil- 
lation parameters (tan 2 #12, Am^) with unprecedented 
precision . Previous analysis of the attainable ac- 

curacy in the determination of the oscillation parameters 



at KamLAND have studied the effect of the time depen- 
dent fuel composition, the knowledge of the flux uncer- 
tainty, the role of the geological neutrinos as well as the 
combined analysis of solar and KamLAND data. These 
studies have been performed in the simplest two-neutrino 
oscillation scheme or equivalently for three-neutrino os- 
cillations assuming a small fixed value of Q 13 

In this letter, we revisit the problem of how the pre- 
cision to which KamLAND should be able to measure 
the solar oscillation parameters, Am 2 2 and #12, is af- 
fected by our ignorance on the exact value of the mixing 
angle #13. At present our most precise information on 
this parameter comes from the negative results from the 
CHOOZ reactor experiment Jl9| , which, when combined 
with the results from the atmospheric neutrino experi- 
ments results into a 3er upper bound sin 2 #13 < 0.06 
|pi| , p2| . To address this question we study how the re- 
constructed range of Am 2 2 and #12 depends on the value 
of the mixing angle #13. We conclude that the recon- 
structed Am 2 2 is very mildly affected by the value of #13, 
while the 6*12 range scales with #13 in a simple way. We 
determine the reconstructed region of solar parameters 
obtained from a given signal once #13 is left free to vary 
below the present bound. We find that the net effect is 
approximately equivalent to that of an uncertainty on the 
overall neutrino flux normalization of up to ~ 10% and 
it should be taken into account once enough statistics is 
accumulated. 

KamLAND is a reactor neutrino experiment located at 
the old Kamiokande site in the Kamioka mine in Japan. 
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It is sensitive to the v e flux from some 10+ reactors which 
are located "nearby." The distances from the different re- 
actors to the experimental site vary from slightly more 
than 80 km to over 800 km, while the majority (roughly 
80%) of the neutrinos travel from 140 km to 215 km. 
KamLAND "sees" the antineutrinos by detecting the to- 
tal energy deposited by recoil positrons, which are pro- 
duced via v e + p — > e + + n. The total visible energy 
corresponds to E e + + m e , where E e + is the total energy 
of the positron and m e the electron mass. The positron 
energy, on the other hand, is related to the incoming 
antineutrino energy E e + = E v — 1.293 MeV up to cor- 
rections related to the recoil momentum of the daughter 
neutron (1.293 MeV is the neutron-proton mass differ- 
ence). KamLAND is expected to measure the visible 
energy with a resolution which is expected to be better 
than a(E)/E = 10%/ VE, for E in MeV 

The antineutrino spectrum which is to be measured at 
KamLAND depends on the power output and fuel com- 
position of each reactor (both change slightly as a func- 
tion of time), and on the cross section for v e +P — > e + +n. 
For the results presented here we will follow the flux and 
the cross section calculations and the statistical proce- 
dure described in Ref. [is{| . We use one "KamLAND- 
year" as the amount of time it takes KamLAND to see 
800 events with visible energy above 1.22 MeV. This is 
roughly what is expected after one year of running (as- 
suming a fiducial volume of 1 kton), if all reactors run 
at (constant) 78% of their maximal power output Jl2| . 
We assume a constant chemical composition for the fuel 
of all reactors (explicitly, 53.8% of 235 U, 32.8% of 239 Pu, 
7.8% of 238 U, and 5.6% of 241 Pu, see ||||). 

The shape of the energy spectrum of the incom- 
ing neutrinos can be derived from a phenomenological 
parametrisation, obtained in p4|], 
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where the coefficients <Xi depend on the parent nucleus. 
The values of a, for the different isotopes we used are 
tabulated in |24|,|l6|]. These expressions are very good 
approximations of the (measured) reactor flux for values 
of E v > 2 MeV. 

The cross section for v e + p — > e + + n has been com- 
puted including corrections related to the recoil momen- 
tum of the daughter neutron in pq | . We used the hydro- 
gen/carbon ratio, r=1.87, from the proposed chemical 
mixture (isoparaffm and pseudocumene) [p^j . It should 
be noted that the energy spectrum of antineutrinos pro- 
duced at nuclear reactors has been measured with good 
accuracy at previous reactor neutrino experiments (see 
fl2"f for references). For this reason, we will first as- 
sume that the expected (unoscillated) antineutrino en- 
ergy spectrum is known precisely. Some of the effects of 
uncertainties in the incoming flux on the determination 



of oscillation parameters have been studied in [ fL6| , and 
are supposedly small. 

In order to simulate events at KamLAND, we need to 
compute the expected energy spectrum for the incoming 
reactor antineutrinos for different values of the neutrino 
oscillation parameters mass-squared differences and mix- 
ing angles. In the framework of three- neutrino mixing 
the v e survival probability depends on the two relevant 
mass differences and the three-mixing angles but it can 
be simplified if we take into account that: 

- Matter effects are completely negligible at KamLAND- 
like baselines. 

- As observed in p), for Am 2 >3x 10 -4 eV 2 , the de- 



termination of Arnr is rather ambiguous. This is due to 
the fact that if Am 2 is too large, the KamLAND energy 
resolution is not sufficiently high to resolve the oscillation 
lengths associated with these values of Am 2 and E v . This 
allows us to consider the higher Am 2 3 and Am 2 3 to be 
averaged. 

Thus, the relevant (energy dependent) electron an- 
tineutrino survival probability at KamLAND is 



P(v e <-» v,,) = sin 4 6<i3 + cos 4 0i 3 



1 - ^ fi sin 2 20i2 sin 2 



(2) 
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where Li is the distance of reactor i to KamLAND in km, 
E v is in GeV and Am 2 2 is in eV 2 , while fi is the fraction 
of the total neutrino flux which comes from reactor i (see 

From Eq. (|3J) we can easily derive the effect of the 
nonzero #13. The energy independent term contains 
the factor sin #13 + cos 4 0i 3 while the energy dependent 
term contains cos 4 0i 3 . Thus the shape of the spectrum 
(this is, the ratio of the energy dependent versus the 
energy independent term) is only modified by a factor 
cos" ^ 1— sin 4 13 . Given the present bound we 

sin 4 #i3+cos 4 #13 10 ^ 

conclude that #13 does not affect significantly the shape 
of the spectrum which is the most relevant information in 
the determination of Am 2 2 . Conversely, the overall spec- 
trum normalization is scaled by cos 4 013-1-2 sin 2 0i3 
and this factor introduces an non-negligible effect. 

In order to quantify the effect of this term we have 
simulated the KamLAND signal corresponding to some 
points in the parameter space (see Table |jj) . Following the 
approach in Ref. |lq ] our simulated data sets are analysed 
via a standard x^unction, 

X 2 (0i 2 ,Am 2 2 ,0i 3 ) = 



E 



Nj (0i 2 , Am 12 , 13 ) ) - T. } : (0i 2 , Am? 2 , i3 ) 



+ V d . . f , 



where Nj(6i2, Am 12 , 013) is the number of simulated 
events in the j-th energy bin which would correspond 
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to the parameters 812 Am 12 813 (see first column in 
Table [| for the values of the 5 simulated points). 
Tj (6 12, Am 2 2 , 0i3 ) is the theoretical prediction for the 
number of events in the j-th energy bin as a function of 
the oscillation parameters. A^in = 12 is the total number 
of bins (binwidth is 0.5 MeV), and the added constant, 
•Wd.o.f j is the number of degrees of freedom. This is in- 
cluded in order to estimate the statistical capabilities of 
an average experiment. An alternative option would be 
not to include the N^.o.t term but to include random sta- 
tistical fluctuations in the simulated data as done in Ref. 
|fj"3"f . We have verified that our results are not quantita- 
tively affected by the choice of simulation procedure. The 
fit is first done for visible energies 1.22 < E V { S < 7.22 
MeV. Note that we assume statistical errors only, and 
do not include background induced events. This seems 
to be a reasonable assumption, given that KamLAND 
is capable of tagging the v e by looking for a delayed 7 
signal due to the absorption of the recoil neutron. There 
still remains, however, the possibility of irreducible back- 
grounds from geological neutrinos in the lower energy 
bins (Evis < 2.6 MeV) To veri fy how this Pos- 

sible background may affect the effect here studied we 
have repeated the analysis discarding the three lower en- 
ergy bins i.e. considering only events with visible energies 
2.72 < E ms < 7.22 MeV. 

We have generated the signal for the five points in pa- 
rameter space listed in Table |[ For the sake of concrete- 
ness we have chosen the five points with different values 
of Am 12 and tan 2 612 distributed within the 3<r allowed 
LMA region from the present analysis of the solar data 
jllj and with 6*13 = 0. For each of the simulated points 
we obtained the reconstructed region of parameters in 
the plane Am 2 2 tan 2 812 by finding the minimum \ 2 and 
then calculating the confidence level (CL) for two degrees 
of freedom assuming three KamLAND- years of simulated 
data. The number of expected events corresponding to 
each of the simulated points with 1.22 < E V i S < 7.22 
MeV is given in Table | and in Table || for 2.72 < E ms < 
7.22 MeV. We have repeated this procedure for the same 
five simulated signals under different assumptions for the 
value of the reconstructed #13. 

In Fig. [|(a) we show the allowed regions in the 
(tan 2 812, Am 2 2 ) plane assuming that we know a priory 
that 9 13 — (which is the simulated value). In other 
words in our minimization procedure we fix #13 = 
in Tj and the only fitted parameters are tan 2 812 and 
Am 2 2 . This case corresponds to the usual two-neutrino 
analysis. Given our \ 2 prescription the best fit recon- 
structed point corresponds exactly with the simulated 
point. The shown regions correspond to la, 2a and 3<r 
for 2 d.o.f (Ax 2 = 2.30, 6.18 and 11.83 respectively). 
Similar regions are obtained if 813 if different from zero 
but assumed to be known so that its value in the simu- 
lated number of events and in the reconstructed one is 



12 



kept to be the same and constant. In the fourth col- 
umn in Table | we list the 3a allowed range for tan 2 
for each of the five simulated points. Table || contains 
the corresponding results from the "conservative" analy- 
sis in which the three first bins have been removed. The 
main effect is a decrease in the statistics which translates 
into slightly larger ranges. We only list the reconstructed 
range in the first octant but, as shown in the figure, due 
to the negligible matter effects, the results from Kam- 
LAND will give us a degeneracy in 812 and an equivalent 
range is obtained in the second octact corresponding to 
tan 2 6*12 — ► 1/ tan 2 8\2- Solar data will be able to select 
the allowed region in the first octant and we show the 3cr 
contours of the LMA region from the latest analysis. As 
discussed in Ref. |18|, if the mixing angle is far enough 



from maximal mixing, the allowed region is clearly sepa- 
rated from the mirror one. 

In order to illustrate the effect of the unknown 6*13 
we have repeated this exercise but now using a different 
value of #13 for the simulated point and the reconstructed 
ones. In Fig. 0(b) we show the reconstructed regions in 
(tan 2 8\2 , Ato 2 2 ) corresponding to the same five gener- 
ated points (which are marked by stars in the figure) but 
using 6*i3 = 12.6° (sin 2 6>i3 = 0.048) in the calculation of 
the expected number of events Tj. From the figure we 
see that best fit reconstructed points (which are marked 
by squares in the figure) as well as the allowed regions 
are shifted in mixing angle with respect to the ones in 
Fig. ^(a) while Amf 2 remains practically unaffected. In 
tables I and H we list the reconstructed ranges of tan 823 
for this academic case. 

The observed shift can be easily understood as follows. 
From Eq. (S), the total number of events is equal for 
different 813 with the condition 



1 - a x sin 2 28 12 = cos 4 13 (1 -ax sin 2 28[, 



(3) 



where a is a number coming from all the detailed in- 
tegration of the oscillating phase factor which depends 
mainly on Am 2 2 and which, for completeness, we also 
list in Table |. For example for the simulated point (1), 

tan 2 !h~2 = 0.37, A^, = 3.7 x 10~ 5 eV 2 and 8^=0, 
the 3(7 reconstructed range, 0.31 < tan 2 0i2 < 0.43, for 
#13 = is shifted using the above relation with a = 0.65 
to 0.25 < tan 2 8 12 < 0.35 for tan 2 8 13 = 0.05 which pre- 
cisely coincides with the values listed in the 5th column 
of table 0. Strictly speaking this scaling is slightly vio- 
lated due to the change of the spectral shape with the 
change from 812 to 8' 12 which also worsens the Xmin f° r 
the reconstructed point. 

Let us finally consider which are the allowed regions in 
the parameter space (tan 2 812, Am|i) taking into account 
that we just know that 813 is below some limit. In order 
to do so one must integrate over 813 in the allowed range, 
or what is equivalent, for each pair (tan 2 812, Am|i) we 
must minimize x 2 (^i2j Ato 21 , tan 2 813) with respect to 
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FIG. 1. (a) Regions of the (Aroli X tan 2 #i2)-parameter 
space allowed by three KamLAND-years of simulated data at 
the la, 2a and 3<r CL, for different input values of Am 21 and 
tan 2 812 (6 13 = 0) and assuming 613 = 0. The stars indicate 
the best fit points (corresponding also to the simulated sig- 
nals), (b) Regions of the (Am|i X tan 0i2)-parameter space 
allowed by three KamLAND-years of simulated data at the 
la, 2a and 3cr CL for $13 = 12.6° (sin 2 13 = 0.048). The 
stars indicate the points used to simulate the signal (with 
fixed 813 = 0) while the squares indicate the reconstructed 
best fit point (with fixed 13 = 12.6°). 
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FIG. 
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la, 
.2 , 



2cr and 3cr CL allowed regions of the 
(Amli x tan 2 6i 2 )-parameter space when 613 is left free to 
vary in the allowed range (0i3 < 13.8° (sin 2 0i3 = 0.057)). 
The regions are obtained for three KamLAND-years of simu- 
lated data and for the same five simulated signals as Fig. 1(a). 
The stars indicate the best fit points (corresponding also to 
the simulated points). 



6*i3 (restricted to be below the present bound). Notice 
that, below the bound, we used a flat probability dis- 
tribution for 6*i3 to keep the analysis just KamLAND- 
dependent. In the future, the combined analysis of solar, 
atmospheric and CHOOZ results with KamLAMD data 
will allow us to include the probability distribution for 
6*13. In Fig. ||we shown the allowed regions for such #13- 
free analysis, where free means allowed to vary below its 
3cr limit, 9 13 < 13.8° (sin 2 6 13 < 0.057)). As seen by 
comparing Fig. || with Fig. [I] the reconstructed regions 
are enlarged and roughly correspond to the overlap of 
the allowed regions for the different fixed values of #13. 
Given our signal generating procedure the best fit re- 
constructed point corresponds to the simulated point in 
(tan 2 612, Am^) and 6*13 = 0. The reconstructed ranges 
can be read from the last column in Tables | and Ta- 
bles ||. Also comparing the results in both tables one 
sees that the presence of this effect is not quantitatively 
affected by not including in the fit the data from the low- 
est energy bins. The main effect being a decrease in the 
statistics which translates into slightly larger ranges. 

We have studied the role of 6 13 in the hypothetical 
case of perfect knowledge of the overall flux normaliza- 
tion and fuel composition. Comparing our results with 
the expected degradation on the parameter determina- 
tion associated with the uncertainty on those two as- 
sumptions 1 16 1, we note that the reconstructed ranges 
obtained in the #i3-free case are very close to those ob- 
tained in the analysis with perfect knowledge of #13 but 
where the overall flux normalization is unknown by about 
~ 10%. As mentioned above, the role of 6*13 is essentially 
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Signal = (tan 2 612, Am^i, #13) 


N ev 


a 


013 = 


013 = 12.6° (0.048) 


0i3 < 13.8° (0.057) 


1 = (0.37,3.7xl0- 5 ,0) 

2 = (0.60,2.5xl0~ 5 ,0) 

3 = (0.70,7.0xl0 _5 ,0) 

4 = (0.25,3.0xl0~ 5 ,0) 

5 = (0.50,2.0xl0~ 4 ,0) 


1022 
802 
1335 
1237 
1321 


0.65 
0.65 
0.56 
0.65 
0.44 


[0.31,0.43] 
[0.48,1.00] 
[0.47,1.00] 
[0.21,0.30] 
[0.37,0.85] 


[0.25,0.35] 
[0.39,0.68] 
[0.35,0.57] 
[0.16,0.24] 
[0.25,0.44] 


[0.24,0.43] 
[0.38,1.00] 
[0.34,1.00] 
[0.15,0.30] 
[0.22,0.85] 



TABLE I. Reconstructed ranges for tan 2 0i2 at 3 a (in the first octant) for different 0i3 (sin 2 0i3 in parenthesis) cases for the 
simulated points listed in the first column. See text for details. 



Signal = (tan 2 0i 2 , Am 2 i, #13) 


N ev 


a 


013 =0 


013 = 12.6° (0.048) 


013 < 13.8° (0.057) 


1 = (0.37,3.7xl0~ 5 ,0) 

2 = (0.60,2.5xl0" 5 ,0) 


756 
695 


0.65 
0.72 


[0.31,0.44] 
[0.41,1.00] 


[0.25,0.36] 
[0.35,1.00] 


[0.24,0.44] 
[0.34,1.00] 


3 = (0.70,7.0xl0~ 5 ,0) 

4 = (0.25,3.0xl0~ 5 ,0) 

5 = (0.50,2.0xl0~ 4 ,0) 


1167 
1004 


0.52 
0.65 


[0.44,1.00] 
[0.20,0.42] 


[0.32,0.59] 
[0.15,0.68] 


[0.31,1.00] 
[0.14,1.00] 


1090 


0.44 


[0.36,1.00] 


[0.24,0.46] 


[0.22,1.00] 



TABLE II. Same as table 1 after removing the lower energy bins from the fit (therefore including only events with 
2.72 < E vls < 7.22 MeV). 



the change of the normalization. This is a larger normal- 
ization error that the 3% expected one in the theoretical 
calculation of the flux from the reactors (induced from 
the /3-spectroscopy experiment at the Goesgen reactor 
p6[). We also find that the uncertainty associated with 
#13 has a larger impact on the determination of the mix- 
ing angle #12 than the expected error associated with the 
fuel composition although, unlike this last one, it does 
not affect the determination of Ani[ 2 . 

Let us point out, that obviously, in order for this effect 
to become relevant it has to be larger than the expected 
statistical uncertainty on the overall flux normalization. 
If we repeat this exercise assuming only one KamLAND- 
years of simulated data, we find very little difference be- 
tween the results corresponding to fixed # 13 = and the 
free-#i3 analysis, as expected since the expected number 
of events would be of the order ~ 400 and the associ- 
ated statistical uncertainty for the overall normalization 
would be comparable with the maximum effect associ- 
ated to #13. 

Summarizing, the KamLAND reactor neutrino exper- 
iment, which is to start taking data very soon, is sensi- 
tive to the LMA region and should be able of measuring 
the solar oscillation parameters (tan 2 #12, Am^) with un- 
precedented precision. In this letter, we have addressed 
the question of the degradation on the determination of 
the oscillation parameters associated with our ignorance 
of the exact value of #13 . We have shown that the deter- 
mination of Ami2 is practically unaffected because the 
effect of 6*i3 on the shape of the spectrum is very small. 
The dominant effect is a shift in the overall flux normal- 



ization which implies that the reconstructed 612 range 
scales with 6*13 in a simple way. As a consequence the 
allowed region of solar parameters obtained from Kam- 
LAND signal will be broader in 6*12. Comparing this 
effect with the ones from the expected uncertainties as- 
sociated with the theoretical error on the overall flux nor- 
malization and the fuel composition, we find that, after 
enough statistics is accumulated, the uncertainty associ- 
ated with #13 may become the dominant source of degra- 
dation in the determination of the mixing angle 0\2 at 
KamLAND. 
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